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SUKHARY 


The  investigation  of  electrochemical  energy  conversion  systems 
has  involved  two  specific  tasks.  The  first  has  dealt  with  measure¬ 
ments  of  the  solubility  of  oxygen  in  organic  liquids  which  might 
provide  the  basis  for  high  energy  batteries.  The  solubility  of 
oxygen  in  y-butyrolactone  is  5.5  ml/100  ml,  in  propylene  carbonate, 

1.7,  in  dimethyl  sulfoxide,  3*7 *  and  in  N-nitrosodimethylardne,  6.6. 

The  physical  method  used  for  these  determinations  gives  results  about 
2.7i°  low  when  compared  with  literature  values  for  the  solubi’ity 
of  oxygen  in  water.  No  change  in  solubility  was  observed  in  butyro- 
laetone  due  to  the  presence  of  added  eloctrolyto  -  up  to  0.8  M 
lithium  perchlorate. 

The  second  task  deals  with  the  mathematical  analysis  of  electro¬ 
chemical  energy  conversion  devices.  Attention  was  given  to  the 
solution  of  kinetic  equations  describing  the  hydrogen-deuterium 
exchange  of  an  hydrocarbon  on  a  catalytic  surface  saturated  with 
deuterium.  Kinetic  expressions  for  a  five  step  successive  reaction 
for  methane  are  solved  by  matrix  techniques.  These  techniques  should 
be  conveniently  handled  by  a  digital  computer  and  the  expressions  can 
be  corrected  for  temperature,  flow  rates,  and  other  experimental 
pa  ranetors. 
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A,  Objective 


In  the  third  progr«*s  report  (l)  on  this  project  the  concept 
snd  theoretic*!  *d  vintage*  of  the  sir  electrode  were  discussed. 

To  Interpret  the  experimental  results  fro*  n ass  transport  con¬ 
siders  tions  it  is  necessary  to  know  the  solubility  of  0^  In 
propylene  carbonate.  To  correlate  cur  results  with  those  obtained 
In  other  systems  It  Is  required  to  know  the  solubility  of  oxygen 
in  organic  solvents  such  ss  y-butyrolaetene,  dimethyl  sulfoxide, 
dine  thy  1  forms  ride,  snd  iJ-nitrosodlrethylanine,  A  survey  of  the 
literature  failed  to  yield  the  pertinent  data  on  the  solubility 
of  0^  in  these  solvents.  For  these  reasons  the  solubility  of  0^ 
was  measured  in  s  rmber  of  organic  solvents. 


£•  Solubility  of  Gases  in  liquids 

There  have  beer,  many  approaches  to  the  explanation  of  gaa 
solubility  in  liquids  including  various  W37S  of  treating  the 
experimental  date.  Sore  of  these  are  discussed  in  the  monograph 
by  Hildebrand  snd  Scott  (2),  None  of  the  approaches  have  been 
universally  applicable. 

One  of  the  r.ore  successful  treatments  of  aolubillty  data 
utilizes  a  plot  of  log  l  (Cstwald  absorption  coefficient*) 
against  r  (solvent  surface  tension)  for  a  given  gas  in  a 


/V  where  v  is  volume  of  gaa  abso; 
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of  solvent  absorbing. 


:d  and  V  is  volume 
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1)  The  solvont  is  degassed.  This  is  accomplished  by  vacuum  distill** 
tlon  using  the  Hestor-Faust  spinning  band  distillation  column.  The  dis¬ 
tillate  is  collected  in  e  flask  which  Is  then  attached  to  the  top  of  the 
absorption  apparatus  via  e  ground  glass  fitting. 

2)  The  U-portion  connecting  the  two  burettes  is  filled  with  pure 
solvont  until  the  liquid  stands  at  the  zero  level  in  both  burettes. 

3)  The  ojygen  gsa  is  flowed  through  the  two  burettes  to  completely  fill 
them  and  to  saturate  the  U-  portion  of  solvent. 

k)  Both  columns  of  gas  sre  returned  to  atmospheric  pressure  and  all 
outlets  to  tho  atmosphere  sre  closed  off. 

5)  The  solvent  in  ths  reservoir  is  slightly  pressurized  with  tydrogen 
to  assume  that  it  will  flow  into  ths  apparatus.  This  is  done  through  s  gas 
inlst  tube  connected  to  the  sample  resorvolr  flask. 

6)  The  solvent  is  dripped  Into  the  absorption  spiral  at  s  rate  of  about 
2.5  ml/minute  and  oollected  in  a  volumetric  flask  as  it  flows  out  of  the 
burette  at  tha  same  rate. 

7)  After  oollsoting  some  appropriate  volume  in  ths  volume  trio  flask, 
such  as  100  ml,  the  burettes  sre  drained  until  the  non-absorption  aids 
reads  0.0  ml.  This  indicates  that  that  side  is  at  atmospheric  pressure. 

Tha  amount  of  solvent  above  0.0  in  tha  absorption  slda  is  tbs  amount  of 

0  absorbed]  this  must  be  sddsd  to  tho  volume  collected. 

2 

K.  Results  and  Disouasion 

On*  important  modification  made  in  the  apparatus  over  that  described 
by  Korrison  and  Biilett  (9)  involved  the  usa  of  a  slight  pressure  head  to 
maintain  a  flow  of  solvent  into  absorbin  burette.  This  pressurizing  gas 
will  be  absorbed  to  sono  degreb  and  then  will  bo  released  into  tha  spiral. 

Thin  obviously  decreases  the  observed  solubility  of  oxygen  in  the  organic 
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liquid  by  some  Magnitude,  If  the  pressurizing  gas  is  oxygen  t  large  negative 
error  would  be  expected.  If  the  pressurizing  gas  wore  N^,  or  a  smaller 
error  would  be  expected  and  the  question  wae  whether  the  error  would  be 
small  enough  to  be  tolerated.  To  investigate  this  effect  and  also  to 
study  the  precision  end  accuracy  of  the  method,  e  number  of  experiments 
were  performed  with  water  ss  the  solvent.  These  values  ere  reported  in 
Table  X.  As  expected,  when  0^  was  used  as  the  pressurizing  ges  the  valuos 
were  quite  low— by  58^.  Although  the  pressurizing  gas  does  not  saturate 
the  water  during  the  time  of  exposure  in  the  reservoir,  considerable  gas  is 
absorbed,  A  series  of  experiments  with  as  the  pressurizing  gss  gave 
results  that  were  7*1?  too  low  while  produced  an  error  of  -  2.l£.  The 
nitrogen  runs  ere  not  replicates  as  certain  improvements  in  technique  were 
Introduced  as  the  experimental  work  progressed.  For  example,  in  the  early 
work  the  elapsed  time  between  degessing  the  solvent  and  running  tha 
axperlment  was  often  several  hours.  Velues  closer  to  the  accepted  litera¬ 
ture  value  were  obtained  when  the  experiment  was  performed  directly  after 
dagassing  through  vacuum  distillation.  Experiments  i«3,  Hi,  and  H2  were  run 
in  this  manner  and  the  average  coefficient,  0.0277,  was  lover  than  the 
accepted  0.028^7  value  (8)  by  2.7£.  The  results  given  below  for  oxygen 
solubilities  in  the  organic  liquids  should  then  be  assumed  to  be  low  by 
about  3/b  assuming  thet  the  solubility  behavior  in  these  orgsnic  solvents  is 
similar  to  that  in  HjO. 

In  Table  II  are  givon  solubilities  of  ojqrgen  in  y-butyrolactone, 
propylone  carbonate,  dimethyl  sulfoxide,  and  K-nitrocodimethylamine. 
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jy  use  of  oth*r  techniques,  theso  expressions  can  be  solved  nore  efficiently. 
b.  lls'-'ir.r.ier. 

Solution  of  the  rate  expressions  for  a  five  stop  successive  reaction 
would  be  desirable  in  order  to  create  a  nodel  for  reactions  such  as  the 
successive  one  site  exchange  of  1/  for  a  in  C.i^.  When  gaseous  CH^  is  passed 
ovar  a  catalytic  surface  which  has  been  exposed  to  such  that  all  available 
■itec  are  assured  to  to  occupied  and  only  cr.o  exchange  of  b  for  H  is  r.ade 

with  each  contact  of  tho  r.oleeulo  with  the  surface,  the  reaction  is: 

t  -'I'v  -  ^  \  '  V.i-»  <♦ 


«  ,  vl  *-  •  v,  \ 


If  it  la  aasuned  that  the  rate  of  increase  of  a  deuterated  species  is  pro¬ 
portional  only  to  tha  nur.be  r  cf  nolo  rules  of  the  lower  deuterated  species 
and  the  probability  of  transition  fron  one  spacie  to  the  next  higher 
deuterated  specie,  first  order  reaction  kinetics  can  be  used  for  each 
transition  bvtveen  specie  and  the  rate  expressions  are: 


(l)  l\V  K.J 


•V\‘  'V'i 


-  \v  i„*  \«.h  -  O'  »• 


Vi-  v.>.i  •  f  . 

\Sll  (  sv»  **>  j'  .*  v'  •  ■'*< 

. J,  . 

.it  '  •  - 

.iV*  '->-1  » 

The  solutions  for  these  expressions  are  given  in  Appendix  A  and  ares 


s  V  .  V..  (■  -  V  1 


‘  ^  Vs*..  V.A-t-.  **<.*<.  i 
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Larson  and  Hall  (12)  assumed  the  expressions  (l)  to  be  applicable 

in  the  H-li  exchange  of  CH^  on  a  deuterated  silica-alumina  surface.  In 

the  absence  of  the  eolutions  to  all  the  rate  expressions,  the  authors 

solved  the  firet  rate  expreseion  for  the  value  of  k  ,  and  assumed  that  the 

ol 

rate  of  exchange  per  hydrogen  atom  is  (?)  k  ,  regardless  of  the  noleculo 

ol 

reacting,  i.e.,  Vo  As.,  /b. .  This  assumption 

impliee  that  the  rate  of  H-D  exchange  ie  constant,  with  a  correction  made 
for  the  nunber  of  H’s  available.  This  assumption  provided  k  values  for  tho 
remainder  of  the  rate  expressions  which  were  then  numerically  evaluated 
from  the  experimental  data.  The  values  calculated  by  tho  authors  using 
thie  technique  were  consistent  with  the  experimental  values. 

Aa  a  verification  of  the  concentration  expressions  (2),  the  concen¬ 
tration  of  each  of  the  deuterated  species  was  calculated  using  the 

author’s  calculated  value  of  k  and  the  assumptions  -V.t  J/\  ^..•W  Vs  A..., 

ol 

-V^-'/nV ,.  The  valculated  values  agree  with  the  experimental  values  within 
etperimental  orrore,  indicating  that  the  assumptions  regarding  the  k 
values  are  reasonable.  However,  with  the  concentration  expressions  (2), 
it  would  not  have  been  necessary  to  make  an}'  assumptions  regarding  tho  k 
values,  ae  each  k  can  be  solved  for  independently  using  these  expressions. 
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Since  tho  probability  of  a  H-D  oxchange  does  appear  to  be 
independent  of  the  specie  reacting,  a  binonlal  distribution  of  the 
deuto rated  species  when  at  equilibrium  is  suggested  where  \>  'M 

J  j*  °h  ?  .  The  use  of  the  binomial  dietributlon 

ia  allowed  ty  inspecting  the  manner  in  which  each  of  the  deuterated 
apecle  is  formed.  For  examplo,  in  ordor  to  form  D^,  there  must  have 
been  two  contacts  of  the  CH  molocule  with  the  surface  in  which  an  H-D 
exchango  occurrod,  and  two  contacts  with  tho  surface  in  which  an  H-D 
exchange  did  not  occur.  If  p=l/4  is  the  probability  of  any  particular 


H  of  CH^  exchanging  on  the  surface,  then  the  probability  of  two 

2  2 

successful  oxchanges  and  two  unsuccessful  exchanges  ia  p  (l-p)  = 

2  2 

(l/4)  (3/4)  .  Since  the  number  of  ways  in  which  four  H  atome  taken  two 


at  a  time  can  be  arranged  is  7,  ^',A:  -  ,  the  total  probability  of 

forming  Dz  is  "v**"  'T>T -  t  ^N\ 

The  oalculated  values  are: 
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Tho  mean  number  of  H  atoms  exchanged  is  4(l/4)  =  1 
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When  these  probability  values  are  expressed  as  a  ratio: 
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The  equilibrium  distribution  of  the  deuterated  species  from  ref (12) 

are: 

C  CW  ^  f.  \ 

O.'.W.jOj  >•« . 

r«. 

•>, .  sc 

0-'°A  ** 

Thus  far  it  has  been  assumed  that  on1;.'  one  H-D  exchange  can  occur, 
when  a  specie  contacts  the  surface.  If  more  than  one  H-D  exchange  can 
occur  on  contact  with  the  surface,  the  following  diagram  represents  the 
ways  in  which  deuterated  species  can  be  formed: 
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.hen  first,  order  reaction  kir.cticr  arc  a  enured t  the  rate  c.xprorrions  ere: 
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The  solutions  to  the  so  expressions  usin-  the  techniques  described  in 
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expressions. 

C.  tlatrix  Solution  of  ldfferenti.al  Venations 
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The  technique  used  to  solve  the  differential  equations  var  that  of 
"conpart-.ent  analysis"  as  introduced  by  iicnuc'r.-ar.n.  (13,1b) .Z ach  of  the  five 
species  of  the  nethar.e  reaction  is  considered  to  be  a  "ce-partnent."  The 
number  of  nolccules  in  each  compartment  is  a  function  of  tine  and  the 
probability  of  a  transition  of  a  nole-culo  into  ar.d  fron  the  cor, partner, t. 

At  the  ncr.cnt  the  nechanisr  cf  transition  between  compartments  in  not 
important,  rather  it  is  desired  to  create  a  mathematical  model  which  "ill 
explain  the  experimental  values.  Th«  model  then  should  dereribe  the 
number  of  nolccules  in  each  compartment  as  a  function  of  the  above  tro 


VAri  abler. 

It  follows  that  this  model  be  ex 
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thr  -  can  l  postulat-d.  If  tV>  c:i,.oriaontnl  data  can  r.ot  be 
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b.  Conclusions 


Exact  expressions  hove  boen  dorivod  for  ll-D  oxchnngo  of  nothane 
on  a  surface  of  adsorbed  deuterium.  While  tho  oxpressicna  are  awkward 
to  handlo,  valid  approximations  are  available  which  will  ease 
calculations. 

From  the  symmetiy  of  the  expressions  for  the  multiple  aito  exchengo, 
it  appears  that  the  expreaeion6  for  other  molcculeo  con  bo  deduced  with 
a  minimum  of  calculations. 

There  are  aevoral  advantages  to  thia  method  of  solving  kinetio 
expressions  over  calculus  methods.  In  particular  the  matrix  operations 
can  bo  conveniently  handled  by  a  digital  computer,  not  only  to  derive 
the  concentration  oxprossiona  but  also  to  substitute  experimental  values 
into  the  rate  expressions.  Also  in  tho  matrix  form,  perturbation 
matrices  can  be  inoluded  to  apply  any  necossary  corrections  for  tempera¬ 
ture,  flow  rates,  etc.,  without  disturbing  the  original  rate  constant 
expressions. 
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